Independent Bond Fluctuation Approximation to the Ground State of Quantum 

Antiferromagnets 
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A simple approach to estimation of the ground state energy of quantum antiferromagnets is 
developed, based on the approximation that quantum fluctuations around different bonds are inde- 
pendent. A transformation of the spin operators to generate bond quantum fluctuations is devised 
and applied to the classical ground state of the S = i Heisenberg model on the square lattice. The 
ground state energy estimates are as good as spin wave theory or slightly better. This simple picture 
of quantum fluctuations might be useful in more complex models. 



The ground state energy of antiferromagnets are low- 
ered relative to the classical energy of the Neel state by 
quantum effects. For the Heisenberg model, spin-wave 
theory!, gives a relatively good quantitative prediction of 
the shift. Here, I show that a simple estimation of the 
quantum contribution at individual bonds gives surpris- 
ingly good predictions of the ground state energy of XXZ 
and XY models. I also devise an explicit transformation 
of the spin operators that generates fluctuations around 
bonds when applied to the classical ground state. The 
energy of the resulting state is calculated on the square 
lattice and found to be close to the simpler estimate. 

The quantum contribution to the ground state energy 
can mathematically be understood as arising from the 
competition between diagonal and off-diagonal terms in 
the Hamiltonian. The lowest energy of the diagonal part 
is found when the system is localized in a particular basis 
state, while minimizing the energy from the off-diagonal 
part requires a superposition of basis states. The Hamil- 
tonian for the XXZ model is 
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(id) 
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where the sum runs over nearest neighbors. Here, I will 
only discuss values of the anisotropy parameter a in the 
range < a < 1. The Heisenberg model is obtained with 
a = 1. If we take as basis the states where each spin i is 
in an eigenstate of sf, the 'zz '-terms in ([1]) are diagonal 
and the 'H — '-terms are off-diagonal. 

I now restrict the further discussion to bipartite lat- 
tices, which can be divided into two sublattices with no 
interactions between spins on the same sublattice. The 
minimum energy of the diagonal terms is for antiferro- 
magnets (J > 0) found with all spins on one sublattice 
having sf = +S ('up'), and all spins on the other point- 
ing the opposite way (s| = — S) ('down'), where S is 
the length of one spin. This is also the classical ground 
state. The energy is — JS 2 per bond. The off-diagonal 
terms all have magnitude aJS. They couple the classi- 
cal ground state to excited basis states in which one 'up' 
spin has been lowered to s? = S — 1 and a neighboring 



the number of nearest neighbors in denoted by Z, each 
of these two spins have diagonal couplings to Z — 1 other 
spins of length S. In the excited basis states, this di- 
agonal energy has increased by 2JS(Z — 1). Even for 
a = 1, this increase in diagonal energy is substantially 
higher than the off-diagonal coupling, so we may expect 
the zero-point fluctuations to be relatively small. I there- 
fore make the approximation that quantum ground state 
fluctuations about different bonds are uncorrelated. To 
calculate the depression of ground state energy per bond 
we only need to take into account the basis state that the 
Hamiltonian of one bond couples to, with the diagonal 
and off-diagonal terms given above. The ground state 
energy per bond Exxz can then be found as the lowest 
eigenenergy of the matrix 



-JS 2 
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aJS —JS 2 + 2JS(Z - 1) 



(2) 



which is 

Exxz = -JS 2 + JS(Z - 1) - JS^(Z-l) 2 +a 2 (3) 
If we write the corresponding eigenvector in the form 



1 

-V 



(4) 



the relative amplitude of the basis state with exchanged 
spins is then 

1 



V = -W{Z-l) 2 + a 2 -(Z-l)) (5) 
a 



The XY model has the Hamiltonian 
Hxy 



(i,3) 



(6) 



It can easily be handled by the method developed here if 
it is considered the 'XZ' model in the chosen basis 2 : 



H XZ = J5>f a j + S ^] 



(7) 
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'down' spin has been incremented to 



-S + 1. If 
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Model 


Reference value 


Spin-wave 


This work 


ID Heisenberg 


-0.443 :p 


-0.4315 


-0.4571 


Square Heisenberg 


-0.3362^ 


-0.329 


-0.3311 


Cubic Heisenberg 


-0.2998.°. 


-0.299 


-0.2995 


ID XY 


-0.318 d 


-0.299 


-0.3090 


Square XY 


-0.2745° 


-0.27 


-0.2707 


Cubic XY 


-0.2640 e 


-0.26 


-0.2625 



a Exacts. 

b Renormalisationl. 

c Self-consistent mean-field 5 . 

d Exacts. 

e Correlated basis functions: 7 . 



TABLE I. Ground state energy (J = 1) per bond of S = § 
antiferromagnetic models calculated by different methods. In 
all cases the classical energy is -0.25. 

The ground state and the diagonal energies are the same 
as for the XXZ model with a = 0.5, but in addition to 

the 'H — '- and ' - + '-terms found in Q, '++'- and ' '- 

terms now appear. However, these new terms destroy the 
state with lowest diagonal energy so they can be disre- 
garded for the calculation of the ground state energy by 
the present method. The remaining off-diagonal terms 
are JS/2, so the estimated ground state energy per bond 
becomes 

E XY = -JS 2 + JS(Z - 1) - JS\J (Z — l) 2 + 1/4 (8) 

It is simple to do the calculation for other anisotropic 
models with antiferromagnetic spin order. Obviously, 
the strongest interactions should be along the z-direction. 
Models with interactions between i.e. S = \ and 5=1 
spins can also be handled by the same general method, 
although new expressions for both diagonal and off- 
diagonal terms must be devised. 

Table Q] lists the ground state energies for the S = \ 
Heisenberg and XY models on selected lattices, with 
J = 1. Shown are estimates of the values from exact 
methods or accurate numerical calculations (see foot- 
notes for references), the results from first order spin 
wave theorji^ and the values found by equation ^ or 
©. These approximations are seen to be as good as the 
spin wave results, in most cases even slightly better (cal- 
culations on the square lattice XXZ model, a < 1, follow 
below). The good agreement suggests that the simple 
picture of uncorrelated bond fluctuations is a reasonable 
representation of local correlations in the true ground 
state. 

The energy estimates §S§ and (J5J do not take into ac- 
count correlations between different bonds. In fact, the 
diagonal energy change 2JS(Z — 1) used in the deriva- 
tion of @ no longer holds when spins in a neighboring 
bond are exchanged. For this reason, one would expect 
the approximation to be best when the amplitude of spin 
exchange is small. Expansion of ([5]) to first order gives 
r\ a 2 /2(Z — 1), and in fact the results in Table U arc 
better the higher the number of nearest neighbours. For 



cases with S > 5, exchange of neighboring spins does 
not flip the neighboring spins completely. Hence, the 
relative change of the diagonal energy is smaller, and 
the approximation can be expected to work better for 
larger values of S. Indeed, for the ID S = 1 Heisen- 
berg model, equation ([3]) with J = 1 gives the ground 
state energy —a/2 ~ —1.414, close to the very accurate 
numerical estimate -1.401-S (the reference gives 9 more 
significant digits). Results in higher dimensions and/or 
for higher values of S are expected to be even better. 

The derivation of ^ consists of taking the classical 
ground state and introducing uncorrelated perturbations 
at each bond. I call the corresponding state the Inde- 
pendent Bond Fluctuation (IBF) state. It is (without 
normalization) 

|IBF) = [][1- ?7S -- S +/25]|0) (9) 

where |0) is the classical ground state and the product 
runs over all pairs of nearest neighbours with index i 
counting 'up'-spins and j counting 'down'-spins. Clearly, 
r\ is the amplitude of the state with spins exchanged at 
one bond, given by ([5]). This is very similar to the Local 
Ansatz approximation developed by Stollhoff and Fulde 9 . 

The relation [sf, s~~] = 2s| holds generally for spin op- 
erators. In the following, the discussion is limited to the 
case S = |, where [sf ,s~~]|0) = |0) if spin i belongs to 
the sublattice pointing up. Therefore, one may treat the 
spin operators as bosonic with s~ as creation operator 
and sf as destruction operator. This is one of the ways 
to express the basic idea behind spin wave theory^. On 
the 'down'-sublattice, sj destroys the classical ground 
state and [sj , sf] |0) = |0). It might be interesting to 
find operators which destroy the IBF state instead, since 
it has an improved description of local quantum correla- 
tions in the XXZ model, sf does the job to zeroth order 
in rj on the 'up'-sublattice, but it can be seen from (|9]) 
that terms of amplitude —77, with a spin pointing up on 
one of the neighboring sites, still remain. These terms 
can be neutralized by adding operators that construct 
the same terms with amplitude +rj: 

4' = 4 + vJ24 ( 10 ) 

where A/i is the set of nearest neighbours to spin i. This 
construction is not ideal, since sf' and sj' on neighboring 
sites don't commute. A more satisfactory operator 
is 

4 w = s f ±2^*^24 ( n ) 

which still destroys the IBF state to first order in 77, but 
also commutes with t ■ on neighboring sites to first 
order in 77. The positive sign in (JTTJ) holds on 'up'-sites, 
the negative sign on 'down'-sites. I now introduce the 
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operator 

L = vJ2( s i4- s tsj) (12) 

where the sum runs over all pairs of nearest neighbors on 
the lattice, and spins of index i are on the 'up'-sublattice 
in |0) and spins of index j are on the 'down'-sublattice. L 
is obviously anti-Hermitian (L< = —L), so the operator 
transformation 

tf = e- L sfe L (13) 
t\ = er L S \e L (14) 

conserves all commutation relations; the t operators ob- 
tained by this transformation are bona fide spin-^ opera- 
tors, albeit delocalized on the lattice. Equation |l3p can 
be expanded 

t+ = S + + [L,,s+] + ^[L,[L,S+]] + ... (15) 
= 4 + ^ ]T s+ + 0{r, 2 ) (16) 

showing that is exactly tf expanded to first order. 
L can alternatively be written 

L = 2ir,J2 ~ = 2^ J2 & x s ^ 

(*>j) (*>j) 

The asymmetric form of L distinguishes the transforma- 
tion developed here from the Local Ansatz method^, and 
is specific for antiferromagnetic interactions. 

The tf operator destroys the |IBF) state to first order, 
and it exactly destroys the state e~ L |0). This state is 
simply the classical ground state in the (slightly) delo- 
calized spin operators defined by (|13M14[) . Its' potential 
as a ground state candidate can be evaluated by calcu- 
lating the energy per bond 

E°(L) = (0\e L H ije - L \0) (18) 

where is any pair of nearest neighbors. I have ob- 
tained an approximate value of E°(L) for XXZ models on 
the square lattice, by expanding e L Hije~ L to high orders 
in a symbolic computation and calculating the value of 
terms with only s z operators (since (0|sf |0) = (0|s?|0) = 
0). In L, terms were included from a part of the lattice 
around the (i, j)-pair large enough to avoid any finite-size 
effects. 

Table HQ lists the results for a = 0.5, 0.8, 0.9 and 1 (in 
units of J). The row labelled 'Reference energy' contains 
low numerical estimates as described in the caption. The 
next row, labelled '£xxz' are calculated using formula 
([3]). The row labelled 'E°(L), nn corr.' gives the results 
of expanding (|18l) to 7th order with the form of L given 
by (|12p and calculating the minimum energy with 77 as a 



variational parameter. The values of 77 are listed below; 
they agree with the values found by (O to within a few 
percent. To assess the convergence of the expansion, it 
can be noted that the biggest contribution to the energy 
from the 7th order term is -0.0002, for a = 1. Hence, the 
accuracy is probably sufficient for comparing the different 
approaches. 

For all values of a, it is seen that E°(L) is a poorer 
estimate of the ground state energy than Exxz, even 
though the calculation of Exxz involves no free param- 
eters. In order to improve the variational result, I have 
tried the calculation with further terms in the definition 
of L. Since the expression (fT2l changes sign under ex- 
change of indices i and j, it is not obvious how to have 
terms with two spins from the same sublattice, as there 
is no natural way to choose the sign for a given pair. 
Therefore, I have only added terms between third-nearest 
neighbours. There are two types of 3rd nearest neigh- 
bours: one type connected by a knight's move (two bonds 
in a row followed by one perpendicular bond), with pa- 
rameter 7721 , the other connected by three bonds in row 
with parameter 7730. Following this notation, the ampli- 
tude of terms involving nearest neighbours is called 7710. 
For the calculation, all terms in (|18p were expanded to 
5th order, and 6th and 7th order terms in 7710 were added. 
The minimum energies obtained in this fashion are listed 
in the row labelled 'E a (L), 3rd nn corr.' and the pa- 
rameters are provided below. The lowering of the energy 
relative to the case with only nearest-neighbour terms is 
modest, just bringing the result in line with Exxz (ex- 
cept for a = 1). Improvements from correlations between 
more distant spin pairs are probably quite negligible, but 
lower energies might be obtained by adding 4-spin terms 
to L. 

Although the energies of states of the type e _L |0) are 
further from the true ground state energies than the sim- 
ple estimate Exxz (not to mention the energies found 
by more elaborate methods), the existence of low- lying 
states of this new, simple form could still be of some 
interest. In particular, the concept of delocalized spin 



a 


0.5 


0.8 


0.9 


1.0 


Reference energy 


-0.2708 


-0.3037 


-0.3183 


-0.3362 




-0.2707 


-0.3024 


-0.3160 


-0.3311 


E°(L), nn corr. 


-0.2706 


-0.3018 


-0.3149 


-0.3292 


V 


0.0837 


0.1329 


0.1481 


0.1629 


E"(L), 3rd nn corr. 


-0.2707 


-0.3024 


-0.3160 


-0.3309 


r/io 


0.0840 


0.1353 


0.1518 


0.1670 


V21 


0.0025 


0.0092 


0.0129 


0.0161 


mo 


0.0007 


0.0037 


0.0056 


0.0073 



TABLE II. Lowest XXZ-model energies and parameters for 
the state e~ L \0) on the square lattice. Reference energies are 
the lowest values from the different methods listed in ref.— , 
except for a = 1, where the number is the same as used in 
Table U 
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operators might be a useful tool to account for quan- 
tum effects on magnetic order in more complex contexts, 
where the most accurate methods are difficult to apply. 
The transformation (flU)) can also be applied to fermions, 
with s + in (JT^J) replaced by the creation operator etc., 



and possibly used to approximate the ground state of the 
Hubbard model or similar models. 

I am grateful to Kim Lefmann for illuminating discus- 
sions and comments on the manuscript. 
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